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Atmospheric Neutrino Fluxes 

Thomas K. Gaisserfl 

'^Bartol Research Institute, University of Delaware, 
Newark, DE 19716, USA 

This talk is a status report on calculations of the flux of atmospheric neutrinos from the sub-GeV range to 
~ PeV. In the lower energy range {E^, < 1 TeV) the primary interest is in using the atmospheric neutrino 
beam to study neutrino oscillations. In the TeV range and above, atmospheric neutrinos are a calibration source 
and background for neutrino telescopes. 



1. INTRODUCTION 

The discovery of neutrino oscillations with at- 
mospheric neutrinos makes it important to know 
the production spectrum of neutrinos as pre- 
cisely as possible in order to infer the proper- 
ties and parameters of the oscillations from the 
data. It also means that the flux of cosmic-ray 
induced neutrinos is much better measured than 
it otherwise might have been. In addition to 
the extensive measurements at Super-K, |ll there 
were important measurements at Soudan |^ and 
MACRO. Q The measurements cover a range of 
energies and techniques, and they see the beam 
from different locations in the geomagnetic field, 
which exposes interesting effects at low energy. 

Because of oscillations, measuring the cosmic- 
ray neutrino beam is an iterative process in which 
the oscillations and fluxes must be understood 
from the same data. Fortunately, calculation of 
the neutrino spectrum at production is straight- 
forward, and it can be checked by comparison to 
measurements of atmospheric muons. Moreover, 
the evidence for oscillations is robust because it is 
based on ratios, which are better known that the 
absolute normalization of the atmospheric neu- 
trino beam. The anomalous ratio of electron- 
like to muon like events reveals a relative deficit 
of at low energy, and the ratio of upward 
to downward multi-GeV events reflects the path- 
length dependence of j/^ oscillations and defines 
a range of Srrr^. M All this is reinforced by the 



low ratio of stopping to throughgoing upward, 
neutrino-induced muons |^,|^ and by the low ra- 
tio of vertically upward to horizontal throughgo- 
ing muons. 1^,^ A consistent pattern of energy 
and pathlcngth dependence emerges that clearly 
points to neutrino oscillations as the explanation. 
Because atmospheric Ve behave normally to the 
precision measured so far, the main effect must lie 
in the v^, <-> Vt sector (or involve sterile neutrinos, 
which are now disfavored). |^,|| 

Here I review the main features of the calcula- 
tion of the atmospheric neutrino beam at produc- 
tion, emphasizing the simple features that pro- 
vide the basis for the evidence for oscillations. 
This talk is based to a large extent on a recent 
review Q. I organize the material here in order 
of increasing energy. 



The atmospheric neutrino flux is a convolution 
of the primary spectrum at the top of the atmo- 
sphere with the yield (Y) of neutrinos per primary 
particle. To reach the atmosphere and interact, 
the primary cosmic rays first have to pass through 
the geomagnetic field. Thus the flux of neutrinos 
of type i can be represented as 



(pp ® Rp < 
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(1) 



'Research supported in part by the U.S. Department of 
Energy under Grant No. DE-FG02-91ER40626 



where is the flux of primary protons (nuclei 

of mass A) outside the influence of the geomag- 
netic field and Rp(A) represents the filtering effect 
of the geomagnetic field. Free and bound nucle- 



2 





M IMIj 

sub 


1 


1 1 Mill 1 
^ \ 




1 nil 




GeV 




^ N 










multi 


/ \ 
/ \ 


through-going 






1 \ 


muons 








GeV 

: \ 1 


1 \ 

1 


\ 

\ 

\ 








/ 




\ 

\ 

\ 












\ 










Stopping 


\ 

\ 








1 i 


\ muons 










\ 


\ 








\ ! l\ 


\ 


\ 










\ 


\ 










\ 


\ 




'J 














I- 




■V-hJ.LLLI ^1^ 


-■+-I.UJ1I h + L- 


1 INI 



3 

■is 



JO"' 10" 10' 10^ 10-' 10' 10" 

Figure 1. Distribution of neutrino energies that 
give rise to several classes of neutrino events. 



ons are treated separately. 

Each of the factors on the right side of Eq. |l| 
is a potential source of uncertainty. Since the 
uncertainties depend on energy, one needs an es- 
timate of the relative importance of different pri- 
mary energies for a given region of neutrino en- 
ergy. Fig. 1 [0 shows the distributions of neu- 
trino energies for four classes of events. Very 
roughly, sub-GeV events, multi-GeV events, up- 
ward stopping muons and upward throughgo- 
ing muons correspond respectively to primary 
cosmic-ray energies of 10^*°-^, lO^-^^"-^, io2 0±o.5 
and 103 "±i GeV. Fig. 2 shows the response func- 
tion in detail for the sub-GeV events, as defined 
at Super-K §. 



2. LOW^ ENERGIES 

The sub-GeV response is at low enough energy 
so that it depends strongly on geomagnetic lo- 
cation and to some extent on the state of solar 
modulation, as illustrated in Fig. 2. The three 
pairs of curves show the distributions of primary 
energy per nucleon that produce sub-GeV events 
under different conditions. These response curves 
are normalized so that the area under each curve 
is proportional to the corresponding signal in the 



absence of oscillations. "S-K downgoing" is for 
events at Super-K induced by cosmic-rays from 
above; "up" is for events at Super-K induced by 
cosmic-rays from below the horizon; and "No- 
cutofF' is the response in the absence of any geo- 
magnetic field effects. The latter situation is re- 
alized in practice for downward events at Soudan 
and SNO, which have low local geomagnetic cut- 
offs. The local geomagnetic cutoff at Super-K is 
high (11.5 GV for vertically downward events), so 
the rate of downward sub- GeV events at Super- 
K in the absence of oscillations would be signifi- 
cantly lower than the rate of upward events. This 
effect is partially washed out by the angle between 
the direction of the charged lepton, which is mea- 
sured, and the neutrino that produced it, and it 
is reversed for muon-like events by the effect of 
neutrino oscillations ||l|. 

Solar modulation may also have a noticeable ef- 
fect on rates in a sufficiently large statistical sam- 
ple. The dashed lines show the maximum effect 
of solar modulation (corresponding to conditions 
during the last two years of data taking at Super- 
K). Neutron monitor records give an appropriate 
measure of solar modulaton because the neutron 
monitor response is similar to the response of at- 
mospheric neutrinos. The shape of the response 
function of a high-latitude neutron monitor is in- 
dicated in Fig. 2. fill 

2.1. East- West effect 

This effect arises from the systematic bending 
of the positively charged primary cosmic rays in 
the geomagnetic field before they reach the atmo- 
sphere to interact. At low geomagnetic latitudes 
the result is a significantly larger geomagnetic 
cutoff for particles arriving from the East than for 
those arriving from the West. Thus the integrated 
cosmic-ray intensity is greater from the West than 
from the East, and this asymmetry is reflected 
by secondary neutrinos and muons. Observa- 
tion of the East- West effect in muons confirmed 
long ago that the primary cosmic rays are posi- 
tive. |l^,|l^ Now observation of this effect with 
neutrinos at Super-Kamiokande fl^ provides an 
important systematic confirmation of the analy- 
sis and interpretation of the evidence for neutrino 
oscillations. 



3 



Single-ring events with momentum in the inter- 
val 0.4 to 3 GeV/c and zenith angles (6) within 
30° of the horizon were compared to calcula- 
tions of the shape of the distribution in az- 
imuth (0). The comparison is independent of os- 
cillation effects to first order because the distribu- 
tion of pathlengths is constant over the 360° band 
of (j) with 60° < 61 < 120°. Thus the agreement 
between observed and expected shapes shows that 
geomagnetic effects are under control in the at- 
mospheric neutrino flux calculations, reinforcing 
confidence in the interpretation of the up-down 
asymmetry of the fiux of muon neutrinos as a 
consequence of oscillations. 



2.2. Primary spectrum 

An estimate |p7|| of the uncertainty in the pri- 
mary spectrum in the light of recent measure- 
ments by BESS in and AMS ^ is ±5% be- 
low 100 GeV/nucleon, which covers the energies 
relevant for sub-GeV events. This estimate is 
based on a fit to BESS and AMS data alone. A 
more proper estimate of the uncertainty in the 
primary spectrum would use all valid measure- 
ments, including those with larger quoted uncer- 
tainties, to estimate the systematic uncertainty in 
the primary spectrum. Excluding an early mea- 
surement which appears to be anomalously 
high, the measurements cover a range of ^ -1-10% 
to —20% below 100 GeV relative to the measure- 
ments of BESS and AMS. 

Figs. 3 and 4 illustrate the uncertainty in the 
neutrino flux from the primary spectrum uncer- 
tainty. The plots show the calculations using 
three fits to the primary spectrum with the same 
interaction model and cascade calculation [pTf . 
Fig. 5 shows the three assumed spectra. The fit 
of Agrawal et al. |^ was made before the mea- 
surements of BESS, AMS and other recent mea- 
surements |p2| , p3| and reflects mainly the earlier 
measurement which has a lower normaliza- 
tion. The two new flts ||l^ are dominated by the 
data from BESS and AMS. The plot shows the 
all-nucleon spectrum, which includes the nucle- 
ons bound in nuclei. 



2.3. Treatment of hadronic interactions 

Atmospheric neutrinos come from decay of pi- 
ons and kaons produced by interactions of cosmic- 
ray hadrons in the atmosphere. The neutrino flux 
depends primarily on the inclusive cross section 
for p -|- air — s- tt^ + X . Although pion produc- 
tion is peaked near xp = because of the steep 
primary spectrum the most important region of 
phase space for the convolution |^ is in the for- 
ward fragmentation region. For pions the range 
0.1 < xp < 0.5 has most weight, while for kaons 
the range extends to somewhat larger xp because 
of p ^ A -I- K+. 

Event generators used in calculations of atmo- 
spheric neutrinos are based on various accelera- 
tor data. Especially important are several spec- 
trometer experiments with protons in the momen- 
tum range 15-30 GeV/c incident on light nuclei. 
Each experiment has statistical and systematic 
uncertainties. In addition, there are systematic 
uncertainties associated with extrapolation into 
unmeasured regions of phase space ||l^. Com- 
parison of calculations made with different inter- 
action models using the same primary spectrum 
gives an indication of the magnitude of uncertain- 
ties that arise from treatment of hadronic interac- 
tions. Such a comparison is shown in Figs. 6 and 
7 made with the primary spectrum of Ref. |2l| 
(heavy solid line in Fig. 5). 

Figs. 6 and 7 from Ref. ||^ display three inde- 
pendent calculations using flve different represen- 
tations of hadronic interactions. The event gen- 
erator Target- 1 was used in the Bartol neutrino 
flux calculation ||2l|. Target-2.1 is a modified ver- 
sion, which is under development for use in three- 
dimensional calculations. The event generator 
has been adjusted to give better agreement with 
the pion data around Xiab = 0.2. Distributions 
of leading nucleons have been adjusted to give a 
better representation of neutron/proton ratios. A 
preliminary description of this work, which is still 
in progress, is given in Ref. HKKM indicates 
the interaction model used in the earlier calcula- 
tion of Honda et al. ; currently |^ they use 
the event generator DpmjetS ||2^. FLUKA refers 
to the three-dimensional calculation of Battistoni 
et al. |2^ . The differences are primarily due to 
the different representations of pion production. 
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rather than to technical differences. 

Differences among the calculations are signifi- 
cantly larger at Soudan than at Kamioka, which 
points to large differences in the treatment of low 
energy interactions (< 10 GeV), which are mostly 
below the geomagnetic cutoff at Kamioka. Forth- 
coming results from the HARP experiment ||2^ 
should be helpful in reducing this ambiguity. This 
should allow a useful, quantitative comparison 
between Soudan and Kamioka, taking account of 
the rather different geomagnetic environments. 



2.4. Three-dimensional calculations 

An approximation common to the calculations 
used to interpret the measurements of atmo- 
spheric neutrinos so far is that the neutrinos fol- 
low the direction of the primaries that produce 
them. This approximation becomes questionable 
at energies low enough so the typical transverse 
momenta of charged pions are comparable to their 
longitudinal momenta, i.e. for Ei, < 1 GeV. 
Another aspect of one-dimensional calculations is 
that the bending of secondary charged particles 
in the geomagnetic field is also neglected. For a 
fully three-dimensional calculation, one needs to 
generate events from an isotropic distribution of 
primaries over the full primary energy spectrum 
(from pion production threshold) at a dense grid 
of locations uniformly distributed over the sur- 
face of the Earth, while accounting for the an- 
gular dependence of the geomagnetic cutoffs at 
each location. In that case only a tiny fraction 
A/R^) of the neutrinos generated will cross 
a detector of area A. Because the radius of the 
Earth i?® ^ vC4, a brute force calculation is a 
challenge. 

So far, two types of approximation have been 
used for three-dimensional calculations. Battis- 
toni et al. neglect bending of secondary par- 
ticles in the geomagnetic field. In this case, a set 
of cascades generated at a single location (with 
no cutoff) can be used to obtain the neutrino flux 
from all directions. This is done by placing the 
detector where each neutrino crosses the depth 
of the detector (once as it enters and once as it 
exits the Earth). The neutrino is then kept if 



the direction of the primary is such that it would 
have passed the geomagnetic field coming from 
the direction required to make the neutrino pass 
through the detector at its true location. The 
other approximation (e.g. ||30| , ^ , ^ ) is to make 
A large. There are also preliminary versions of 
more ambitious calculations ||33| , ^ . 

An excess of low energy neutrinos from near 
the horizon is a characteristic feature of three- 
dimensional calculations. |^ The origin of this 
geometrical effect is discussed by Lipari. Q Be- 
cause it occurs only for low energy neutrinos, 
however, the excess is washed out by the broad 
angular distribution of the the charged lepton rel- 
ative to the neutrino direction. 

Associated with the horizontal excess is a small 
change in the pathlength distribution of low en- 
ergy neutrinos. Quantitatively, the average path- 
length of neutrinos arriving from just above the 
horizon (80° < 9 < 90°) is about 10% lower in 
the 3D calculation than in the ID approxima- 
tion for El, — 0.3 GeV. The suppression for 
Ey = \ GeV is < 5%. This effect, which has not 
yet been included in analysis of Super-K, could 
be expected to increase the inferred 5m? slightly 
(of order ^ 1% to keep the factor 5m? x LjEy 
constant while decreasing L by ^ 10% for ^ 10% 
of sub-GeV events). 

During its test flight on the Space Shuttle, 
AMS mapped the energy spectrum of low-energy 
primary cosmic rays around the globe. The data 
show a substantial flux of protons with energies 
below the local geomagnetic cutoff. |3^ These 
are secondary protons produced by interactions 
of cosmic-rays above the cutoff that enter the 
atmosphere at large zenith angles. Depending 
on their orientation relative to the geomagnetic 
field, a large fraction of the sub-cutoff secondaries 
curve back into space. Such albedo particles 
may remain trapped for several cycles before they 
re-enter the atmosphere. As a consequence, the 
rate per area-solid angle at which they re-enter is 
correspondingly lower than their flux at 380 km. 
These features have been evaluated quantitatively 
by Zuccon el al. |^ . Since both the flux of re- 
entering subcutoff particles and their pion multi- 
plicity is lower than for primaries above the cut- 
off, their contribution to the neutrino flux should 
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be small. It is interesting to note that all sub- 
cutoff particles are included in the ID calculations 
because all secondaries are assumed to follow the 
direction of the incident primary. In this respect, 
the ID calculation overestimates slightly the con- 
tribution of sub-cutoff particles since a few may 
interact on the way up rather than upon re-entry. 

Another consequence of bending of secondary 
charged particles in the geomagnetic field, in this 
case involving muons, is a second-order system- 
atic effect on the calculation of the East- West 
effect. The asymmetry is increased for and 
reduced slightly for f^. Q 

3. HIGHER ENERGY 

For Ev> \ GeV the energies of the parent cos- 
mic rays are sufficiently high that geomagnetic ef- 
fects and solar modulation become unimportant. 
In addition, the direction of the charged lepton is 
better aligned with the direction of the neutrino. 
Thus for multi-GeV events and neutrino-induced 
muons, in the absence of oscillations, one sees the 
characteristic enhancement of the neutrino flux 
from near the horizontal. This is related to the 
"secant theta effect" familiar from atmospheric 
muons. It is a consequence of the enhanced de- 
cay probability for pions at large zenith angles. 
The corresponding effect for neutrinos is nicely 
illustrated in the plots of angular distributions of 
neutrino events from Super-K. j^] 

3.1. Angular dependence 

An analytic approximation for the differential 
flux of -\- j/p from decay of pions and kaons dis- 
plays some important features the neutrino flux. 



where 



{I-ZnnKi + I) 



1 + -Bjriy COS OE^/eT, 



0.635 



dN 
dE^ 



1 + Bku cosOE^/tK 



(2) 



A-kE, 



-(7+1) 




(3) 



is the differential primary spectrum of nuclcons 
of energy Sq. The neutrino flux in Eq. ^ is 
proportional to the primary spectrum evaluated 



at the energy of the neutrino. The constants 
Ti = nr^/m'^ for i ~ {tt,K), and the constants 
Bi depend on hadron attenuation lengths as well 
as decay kinematics [p9| . 

The critical energy for pions is e^r ~ 115 GeV, 
while for kaons ck « 850 GeV. For E^ < 
e7r/cos6' the neutrino spectral index is the same 
as that of the primary spectrum. For Ei, > 
e^/cos0 the spectrum steepens, first for neutri- 
nos from pion decay and at higher energy [E^, > 
eK / cos 9) for neutrinos from decay of kaons. As 
energy increases so does the horizontal enhance- 
ment from the zenith angle dependence in the de- 
nominator of Eq. Ij. Fig. 8 shows the fraction 
of atmospheric neutrinos and atmospheric muons 
that come from pions and from kaons at two an- 
gles. Above ~ 100 GeV, kaons dominate as the 
source of neutrinos. It is the larger value of £k^ 
together with the form of the kinematical factor 
in the numerator, that determines the behavior 
shown in Fig. 8. 



3.2. Flavor and charge ratios 

The primary source of electron neutrinos is 
muon decay, which also gives an equal contribu- 
tion of muon neutrinos. For i?^ > 2 GeV/ cos 9 
muons on average reach the ground before they 
decay. For 9 > 70° the curvature of the Earth 
limits the pathlength in the atmosphere to a few 
hundred km and muons with i?^ > 30 GeV typ- 
ically reach the ground before decaying. For E^, 
more than ~ 100 GeV/cos6', the contribution 
from muon decay practically vanishes and the ra- 
tio of Vf^jv^ approaches a small value with 
most of the v^. coming from i^'es decays. Some- 
where above a TeV, charm decay dominates Vf^. 

Atmospheric neutrinos constitute the primary 
background, as well as a useful calibration source, 
for neutrino telescopes. The characteristic angu- 
lar and energy dependences of the fluxes of and 

described here should be useful for calibration. 
As an example, Fig. 9 shows the ratio of hori- 
zontal to vertical neutrinos for two flavors from 
Ref. plf The angular asymmetry of muon neutri- 
nos increases steadily with energy to a ratio ^ 2 
in the TeV range. The large horizontal to vcrti- 
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cal ratio for Vg, and its energy-dependence above 
100 GeV should be a distinctive calibration sig- 
nature in detectors large enough to measure the 
low rate of i/e-induced cascades. 

The ratio v^/v^ is also of interest. As mag- 
netized detectors such as MINOS |40 begin to 
operate, it will be possible to measure the charge 
ratio of atmospheric neutrino- induced muons. At 
the high energies the kaon contribution will be 
of great importance for evaluating the u/v ra- 
tio, which, together with the cross sections (Ti, > 
CTp, determines the muon charge ratio. Because 
the contribution of is larger than that of 
< 1. 



3.3. Prompt neutrinos from charm 

Leptonic decays of short-lived charmed hadrons 
gives rise to a third term (not shown) on the 
right side of Eq. |[ of the same form as the 
term for the kaon contribution. The spectrum- 
weighted moment .^Ar^charm is much smaller, but 
echaimXLlO^ GeV, so this component dominates 
the neutrino spectrum at sufficiently high energy. 
There are significant uncertainties in inclusive 
cross sections for hadronic production of charm 
and hence for the prompt neutrino flux. Be- 
cause of the steep cosmic-ray energy spectrum, 
charm production in the forward fragmentation 
region is likely to make an important contribution 
to the flux of prompt leptons, in addition to pro- 
duction via central QCD processes. Fig. 10 shows 
a global view of the spectrum with an estimate of 
charm production from Ref. [4^ . In this quark- 
parton recombination model, charm becomes the 
dominant source of i/g for E^>2 TcV and of j/^ for 
E^>100 TeV. Prompt neutrinos have an isotropic 
angular distribution. The increasing importance 
of prompt h'e contributes to the rapid decrease of 
the ratio of horizontal to vertical electron neutri- 
nos above 100 GeV in Fig. 9. 
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E (GeV) 

Figure 2. Response functions for sub-GeV neu- 
trinos under several conditions (see text). The 
three pairs of curves show response functions for 
the distribution of neutrino energies in the sohd 
curve (S-K i?^), which is the same as the left- most 
curve of Fig. 1. 
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Figure 4. Comparison of neutrino flux for three 
primary spectra for the location of Soudan. 
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Figure 3. Comparison of neutrino flux for three 
primary spectra for the location of Kamioka. 



> 



T 10" 



in 

T3 



10" 



10" 



All nucleon flux 



Agrawal at al. 1996 
New flux-1 
New flux-2 




10' 



-kin 



10-^ 
(GeV) 



10'' 



Figure 5. Three fits to the spectrum of cosmic-ray 
nucleons. 
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Figure 6. Comparison of atmospheric neutrino 
flux calculations for the location of Kamioka. 
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Figure 7. Comparison of atmospheric neutrino 
flux calculations for the location of Soudan. 
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Figure 9. Ratio of horizontal (0. < |cos^| < 
0.375) to vertical (0.675 < |cos^| < 1.) atmo- 
spheric neutrinos. 
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(GeV) 

Figure 10. Global view of the atmospheric neu- 
trino spectrum integrated over all directions. The 
upper solid line shows the total flux of f ^ + Pf^ in- 
cluding the effects of oscillations. Above 100 TcV 
this flux is dominated by prompt neutrinos. Bro- 
ken lines are Ve + De, with the crossover to prompt 
around 1 TeV. Shaded area indicates solar neu- 
trinos. 



